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Geometric frustration is an approach to the glass transition based upon the consideration of locally favoured 
structures (LFS), which are geometric motifs which minimise the local free energy. Geometric frustration 
proposes that a transition to a crystalline state is frustrated because these LFS do not tile space. However, 
this concept is based on icosahedra which are not always the LFS for a given system. The LFS of the popular 
Kob-Andersen (KA) model glassformer is the bicapped square antiprism, which does tile space. Such an 
LFS-crystal is indeed realised in the AI 2 CU structure, which is predicted to be a low energy state for the 
KA model with a 2:1 composition. We therefore hypothesise that upon changing the composition in the KA 
model towards 2:1, geometric frustration may be progressively relieved, leading to larger and larger domains 
of LFS which would ultimately correspond to the AI 2 GU crystal. Remarkably, rather than an increase, upon 
changing composition we find a small decrease in the LFS population, and the system remains impervious 
to nucleation of LFS crystals. We suggest that this may be related to the composition of the LFS, as only 
a limited subset are compatible with the crystal. We further demonstrate that the AI 2 CU crystal will grow 
from a seed in the KA model with 2:1 composition and identify the melting temperature to be 0.447(2). 


I. INTRODUCTION 

Understanding the glass transition remains one of 
the major outstanding questions of condensed matter 
physicsii^. Broadly speaking, theories can be classified 
into those which propose a transition to an “ideal glass” 
state at finite temperature^^— and those which propose 
that there is no thermodynamic transition but that struc¬ 
tural relaxation becomes progressively slower upon cool¬ 
ing towards absolute zero^. Since at some point a glass 
forming liquid cannot be equilibrated on an experimental 
timescale, it is hard to establish whether its structural re¬ 
laxation time Ta would diverge at finite temperature (in 
the case of an “ideal glass”) or whether it it would re¬ 
main finite (although large) until temperature drops to 
absolute zero. Fits to experimental data are insufficient 
to discriminate either scenario conclusively^. 

Although the absolute nature of the glass transition 
remains unclear, in recent years, structural change has 
been measured in a number of materials approaching dy¬ 
namical arrest— S. Structural change is necessary in the 
case that there is a thermodynamic transition, but it is 
not excluded by the no-transition scenario. This struc¬ 
tural change is hard to detect in two-point correlations 
such as the radial distribution function g(r) but can be 
detected using higher-order contributions such as three- 
body correlationsiiii^. Other methods include identify¬ 
ing geometric motifs which minimise the local free en¬ 
ergy, called locally favoured structuresiSiHiii and order- 
agnostic methods where the nature of the structural cor¬ 
relations are not specifiedi^^—. 

The theory of geometric frustration makes explicit ref¬ 
erence to locally favoured structures (LFS)^. In d=3 
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the archetypal LFS is the icosahedron identified by 
Frank when considering the monatomic Lennard-Jones 
modeli^. Icosahedra do not tile Euclidean space, thus 
inducing frustration and preventing crystallisation. The 
situation is rather different in curved space where 600 
perfect (strain-free) tetrahedra comprising 120 particles 
can be embedded on the surface of a four-dimensional hy- 
persphere2i“— . Each particle in this 4d Platonic solid or 
“polytope” is at the centre of a 12-particle icosahedrally 
coordinated shell. Mapping from curved space back to 
Euclidean space, introduces defects in the tiling of icosa¬ 
hedra which leads to frustration and the avoidance of a 
transition to a fully icosahedral state. Thus the degree 
of deviation in curvature from that of the fully icosahe¬ 
dral state corresponds to the strength of frustration. For 
small deviations, frustration is weak and large domains of 
icosahedra can form. Further deviation in the curvature 
corresponds to strong frustration and smaller domains of 
icosahedra. 

Locally favoured structures other than icosahedra are 
also possible. In d=2, the LFS is often a hexagon which 
does tessellate in 2d space. This means that curving 
space induces frustration rather than relieving it as in 
d=321“— . It has been shown that LFS domain sizes can 
be very large in the case of weak frustration; further¬ 
more, the degree of curvature controls the domain size as 
predicted by geometric frustration theor y^^i^^ . Again, in 
d=3 the LFS need not be the icosahedroniS but may be 
other motifs depending which system is considered22ri2i. 
A number of investigations have been made into LFS 
domain size in d=3, with mostj^^— (but not all^“— ) 
finding that the characteristic lengthscale of the domains 
is small. This implies that in 3-dimensional Euclidean 
space, geometric frustration is strong^!. 

The idea of curving space to relieve geometric frustra¬ 
tion is predicated upon the LFS being icosahedral^i^^— . 
In general this is not the only possible method of relieving 
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frustration. Examples include the use of polydispersity to 
control frustration in systems where the locally favoured 
structure was crystalline^i*^ though others found that 
the LFS in weakly polydisperse hard spheres was based 
on five-membered rings with strong geometric frustra- 
tio n^'^ 1 ^^ 1 ^^ . A number of glassformers have been found 
to crystallise, often into complex crystal structures which 
incorporate the LFS — . Demixing into simpler crystals 
such as FCC is also possible^, as is the formation of 
substitutional FCC lattices^ 

In this work we consider the Kob-Andersen (KA) bi¬ 
nary Lennard-Jones model, a popular 3d model glass- 
former in which the LFS is a bicapped square antiprism 
of 11 particles [see Fig. [T] (a)p2i>^. We believe that the 
structure of the LFS is related to the non-additive in¬ 
teractions between the two species in the KA model. 
Supporting this hypothesis, certain metallic glasses with 
non-additive interactions also exhibit considerable popu¬ 
lations of bicapped square antiprism o^°i'^^'^^ . The struc¬ 
ture of the LFS in the KA model is interesting in the 
context of geometric frustration because it tessellates in 
Euclidean space into a crystal of the form Al 2 Cu^. Al¬ 
though the full phase diagram for the KA model has yet 
to be determined, some state points have been consid¬ 
ered. The 1:1 mixture crystallises readily to a CsCl struc¬ 
ture^!^ while the 4:1 composition is predicted to phase 
separate into coexisting face centred cubic lattices^. Ev¬ 
idence in support of these scenarios has been found upon 
quenching, but around the 2:1 composition no evidence 
of either FCC or CsCl was found^. Here we primarily 
consider the 2:1 composition and in particular the role of 
the AI 2 CU structure. Previous work with this composi¬ 
tion found some change in local geometry in the super¬ 
cooled liquid — but focused on smaller structures than 
the bicapped square antiprism which was later identified 
as the LFS^i^s. 

The AI 2 CU crystal offers a new means to control geo¬ 
metric frustration in the Kob-Andersen model. Rather 
than curvature, the composition of the mixture can be 
modified between the 4:1 of the normal mixture to 2:1, 
and the size of the LFS domains can be measured to 
determine the degree of frustration. We begin by consid¬ 
ering the normal 4:1 mixture, measuring the LFS pop¬ 
ulation and investigating the domain size. Comparing 
this to the 3:1 and 2:1 compositions for a given degree 
of supercooling, we unexpectedly find a small decrease in 
the LFS population as the composition is changed, rather 
than an increase. We will show that this appears to be re¬ 
lated to the need to consider the composition of the LFS 
rather than just the number of LFS in the system. In 
moving from the 4:1 composition to the 2:1 composition, 
we find a significant increase in the population of LFS 
which are compatible with the AI 2 CU crystal, but this 
remains a fraction of the total LFS population. We find 
no evidence for nucleation of the AI 2 CU crystal but we 
are able to grow it from a seed and estimate the melting 
temperature as 0.447(2). 

This paper is organised as follows. After describing 



FIG. 1 . (a) The geometric form of a bicapped square an¬ 

tiprism, the LFS for the Kob-Andersen binary mixture, (b- 
c-d) The three most common configurations of the LFS in a 
KA liquid with the larger A particles in white and the smaller 
B particles in purple; (b) AiqBi, (c)A 9B2 and (d)A8B3. 
Only the AgBs form tessellates to form a bulk ALCu crys¬ 
tal. Centre-to-centre distances between the particles defining 
the principal rotational axis of energy-minimised clusters are 
reported next to each bicapped square antiprism’s type. 


our simulation and analysis methodology in section|TTl we 
present details of the dynamical behaviour of our simula¬ 
tions in section ITlI A1 We then show that the composition 
of the LFS is altered by modifying the composition of the 
KA model in section llll B1 In section llll Cl we proceed to 
show that the 2:1 composition can undergo seeded crys¬ 
tal growth in MD simulations to give a crystal of the 
form AI 2 CU. We discuss and interpret our findings in 
section m and conclude in section El 


II. METHODS 

A. Model 

The classical Kob Andersen (KA) binary Lennard- 
Jones mixture is a model glassformer composed of 80% 
large particles (A) and 20% smaller particles (B) and is 
designed to be homologous to the metal alloy, Ni 8 oP 20 j 
a metallic glass^^*^. The Lennard-Jones parameters are, 
o'AA = O', aAB = 0.8a, asB = 0.88cr, eaa = £, 
Sab = l-5e and ebb = O.Se. The particles have the 
same mass, tua = rns- The interactions are truncated 
at = 2.5a, = 2.0a and = 2.2cr. 

In the present work, we consider the usual 4:1 compo¬ 
sition as well as two additional mixtures with the same 
interactions but modified ratios of A to B particles (3:1 
and 2:1). The units used throughout this work are re¬ 
duced Lennard-Jones units, with respect to the A parti¬ 
cles. This means that length is measured in units of a, 
energy (E) in units of e, density (p) in units of Na^/V, 
time (t) in units of ^Jma'^ jE, pressure (P) in units of E/a^ 
and the Boltzmann constant (fc{,) is set to 1. Through¬ 
out we used the LAMMPS molecular dynamics package^ 
and set pressure P = 0. 
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B. Structural Analysis 

In order to detect locally favoured structures in the 
considered systems, we employ the Topological Cluster 
Classification (TCC) algorithm. The TCC is described 
in detail elsewhere^ but in summary, it identifies a neigh¬ 
bour network for all of the particles in the system using 
a Voronoi tessellation and then determines shortest path 
relations for 3, 4 and 5 membered rings. From these ba¬ 
sic structures, more complex structures are identified in a 
hierarchical fashion. The structure that we are interested 
in identifying is the LFS, the bicapped square antiprism. 
In the language of the Voronoi Face Analysis^ this corre¬ 
sponds to the (0,2,8) class^. Because we are considering 
a binary mixture, the composition of the LFS may be 
significant. The TCC is able to detect the composition 
of structures that it recognises in terms of the number of 
“A” and “B” particles that make up the cluster. We find 
that there are three main LFS compositions. All have a 
small B particle in the centre and differ in whether the 
other two positions along the principal rotational axis (or 
spindle) of the structure are occupied by A or B particles 
(see Fig. [1]). There is some small variation in structure, 
with a B particle occasionally occupying a non-spindle 
position, though the three structures shown account for 
95% of the LFS in a 2:1 KA liquid at T = 0.44. We refer 
to these three structures as the AiqBi, the A 9 B 2 and the 
AgBa clusters. 


III. RESULTS 


A. Slow Dynamics 

In order to characterise the relevant time scales of the 
KA mixtures, we have determined their relaxation times 
over a range of temperatures (see Fig. [2^.). We per¬ 
formed simulations oi N = 3375 particles at P = 0 in 
the isobaric-isothermal ensemble for three different mix¬ 
tures: the 2:1, 3:1 and 4:1 compositions. After a quench 
from a random (infinite temperature) configuration the 
systems were equilibrated at each temperature and from 
these equilibrium simulations, the intermediate scatter¬ 
ing functions were computed, providing an estimate for 
the a-relaxation times r^. For each temperature, the 
equilibration time is at least IOOtq,. 

The non-Arrhenius growth of the relaxation times as a 
function of the inverse temperature can be fitted with the 
conventional Vogel-Fulcher-Tamman (VFT) equation 


r„(T) = Too exp 


DTo 

AT-To)_ 


( 1 ) 


which allows for the computation of the fragility D and 
the ideal-glass transition temperature Tq, at which the 
relaxation times appear to diverge. These values are 
reported in Table U and show that the change in com¬ 
position leads to a significant increase in the transition 


temperature Tq when moving from the 4:1 to the 2:1 com¬ 
position accompanied by an increase in fragility. With¬ 
out entering into debate about the physical significance 
of the VFT fiti, we use Tq to rescale our data for conve¬ 
nient comparison of the degree of supercooling between 
the different compositions. As discussed in the following 
section, this change in composition is related to impor¬ 
tant modifications in the structure of the liquids. 


VFT 4:1 3:1 2:1 

T7 0278 0 300 0336 

Too 0.122 0.132 0.142 

D 2.81 2.46 2.10 


TABLE I. Parameters of VFT fits to relaxation times for KA 
liquids of different compositions. 


B. Locally Favoured Structures 


In simulations of the KA mixture, the population of 
LFS is inversely related to temperatur o^^i^^ . Fig. [JDb) 
shows the proportion of particles in a KA liquid that 
are in a cluster of LFS at P = 0 for the three mixtures 
considered in this study. 

We observe that the three mixtures have similar LFS 
concentrations and that, upon rescaling the temperature 
with respect to the fitted ideal glass transition temper¬ 
ature Tq, the 2:1 mixture shows fewer structures for a 
given degree of supercooling [see inset of Fig. [2l[b)]. The 
expected trend is that modifying the system such that 
the LFS is compatible with the crystal should increase 
the population of bicapped square antiprisms but our re¬ 
sults are contrary to this. 

We can consider the relation between structure and dy¬ 
namics by observing the relation between the LFS con¬ 
centration and the relaxation time of the system, see 
Fig. [21:. In this case, it is possible to observe the slowing 
of the dynamics induced by the change in the stoichiome¬ 
try. For an LFS concentration of approximately 10% the 
relaxation time of the 2:1 mixture exceeds that of the 4:1 
mixture by an order of magnitude. Therefore, it appears 
that the slowing dynamics of the system is correlated to 
a change in the composition of the locally favoured struc¬ 
tures, and in particular to an increase in the proportion 
of AgBa clusters. 

Figure Ela) shows the variation in LFS structure for 
the three different KA mixtures that we have considered 
at a temperature T = l.dlTo- Next to it, in Fig. [SIj, we 
report the ground state potential energies of the same 
clusters, evaluated in isolated configurations using an ef¬ 
ficient basin-hopping method^i^. The composition that 
best minimizes energy in the bicapped square antprism is 
the AgBg cluster. We note that it is this cluster which is 
compatible with the AICI 2 crystal, the other two are not. 
Despite its low energy, the AgBg cluster is almost absent 
in the classical 4:1 mixture in fact, which is dominated 










4 



FIG. 2. (colour online) (a) Angell plot for KA liquids with different ratios of large to small particles at P = 0. Black lines 
indicate fits according to Eq. [T] (b) The fraction of particles that are in a cluster of LFS for different ratios of large to small 
particles as a function of the degree of supercooling Tq/T where Tq is reported in Table |T] In the inset, the same fractions 
plotted against the reduced temperature, with no rescaling. The relaxation times and the fraction of LFS are compared in 
panel (c), where an important slowing down is observed for the 2:1 composition. 


by the highest energy cluster AiqBi. Increasing the pro¬ 
portion of B particles in the mixture through the 3:1 to 
the 2:1 mixture leads to a signihcant increase in the rel¬ 
ative amount of AgBa to (« 25%) of the LFS. Although 
these energetically optimal AgBg clusters do tessellate, as 
demonstrated in Fig. SI they are never clustered together 
in the fluid mixtures and during our simulations never at¬ 
tain a critical size sufficient to form a crystalline nucleus. 
This aside, they do optimize the local energy landscape 
and their presence in the 2:1 mixture correlates with the 
exceptional slowing down observed in Fig. [^. 


■A-ioBi ■AgB2 \AgB3 



3 Ratio A:B b 

FIG. 3. (colour online) (a) Composition of LFS clusters in 
KA mixtures of varying A:B particle ratios at T = IAITq. (b) 
Energies of isolated LFS clusters of different compositions, in 
reduced units. 


C. Crystallisation 

Given that we see no crystalline order or even signifi¬ 
cant increases in the size of connected LFS domains su¬ 
percooled liquids, we enquire as to the stability of the 
AI 2 CU crystal. Motivated by the work of Fernandez and 


Harrowell^ who considered the crystallisation of a CsCl 
crystal from a KA mixture using a crystal seed, we con¬ 
ducted seeded crystal growth molecular dynamics simu¬ 
lations in the isobaric-isothermal ensemble at P = 0 with 
a 2:1 ratio KA mixture. 

First we determined the density of the bulk crystal 
in the A^Cu conformation at zero pressure and constant 
temperature for all the low temperatures studied in a sys¬ 
tem of 3375 particles. Then, simulations were performed 
by taking a crystal of A^Cu consisting of 20736 particles 
at the previously determined density and fixing a 1500 
particle cubic region in the centre of the simulation box. 
The surrounding particles were equilibrated at T = 2 and 
then quenched to low temperatures. In this way we ap¬ 
proximate a 2:1 KA liquid surrounding an ideal A^Cu 
seed. The central seed was then unconstrained and the 
growth process tracked in the isothermal-isobaric ensem¬ 
ble for a time equivalent to 500rQ, , where the relaxation 
time refers to the corresponding 2:1 KA liquid at P = 0 
(see section ITlI All . 

To determine the size of the crystalline region, we 
coarse-grain by considering only the centres of the bi- 
capped square antiprisms. Using this coarse-graining, 
we define domains of LFS using an agglomerative single 
linkage hierarchical clustering method^, identifying the 
largest domain as the growing crystal. 

In this way, we detect the stable crystalline structure 
from the LFS belonging to the liquid region, as shown 
in the two snapshots in Fig. [S] From direct inspection 
of the configurations, we notice that the crystal grows 
in a cylindrical fashion, with the main axis aligned with 
the primary rotational symmetry axis of the bicapped 
square antiprisms. Few new clusters are added to the 
lateral surface of the cylinder throughout the simulations, 
suggesting a strong anisotropy in the surface tensions of 
the different crystalline orientations. 

We considered this system for a range of temperatures 













5 



FIG. 4. (colour online) Network of LFS for (a) 4:1, (b) 3:1 and (c) 2:1 Kob-Andersen liquids at T = 1.41To. Purple particles 
indicate clusters of stoichiometry AioBi, green are A 9 B 2 and white are AsBs. 



FIG. 5. (colour online) Snapshots of the system at differ¬ 
ent times: (a) T = 0.425, immediately post quench and (b) 
T = 0.425, after 500ra. Centres of crystalline LFS are repre¬ 
sented in red, centres of non-crystalline LFS are represented 
by crosses. 



0 -^^^ 

0 100 200 300 400 500 


t/Ta 


T = 0.42 — 0.46 which spanned from crystal growth to the 
dissolution of the seed into the melt. The simulation was 
performed 5 times at each temperature with randomised 
initial particle velocities, the average crystal growth is 
plotted in Fig. [6] From this series of non-equilibrium 
simulations we infer the temperature Tmeit at which the 
growth rate is zero. 

For temperatures T < 0.45 we observe the melting 
of the initially seeded structure (see Fig. [SI). This al¬ 
lows us to estimate positive and negative growth rates 
via linear fits to the growth curves (Fig. |7|) and from 
this determine the melting temperature of the crystal as 
r^eit = 0.447(2). 

The composition of the crystal and the liquid was mea¬ 
sured in terms of A and B species at T = 0.42. The crys¬ 
tal has a ratio of large to small particles of 2.35 which is 
close to the value of the 2:1 stoichiometry of the crystal, 
the over-representation of A particles is likely due to the 
the composition of the particles at the boundary between 
the crystal and the liquid being higher in A particles. 

The LFS clusters in the liquid however have a ratio of 
4.15 large particles to each small. This is close to the 4.23 
expected for a 2:1 KA liquid, the composition of which is 
shown in Fig. [3] This suggests that the structure of the 
LFS in the liquid is not strongly affected by the presence 
of the crystal throughout the simulation. Similarly, the 


FIG. 6 . Time evolution of the number of bicapped square 
antiprism centres in the largest crystalline cluster for different 
temperatures. Time t — 0 corresponds to the instant at which 
the initially frozen crystalline nucleus is freed and begins to 
thermalize. Each temperature is an average of 5 simulations. 

percentage of particles in LFS clusters in the liquid re¬ 
gion is about 10.5% which is broadly in line with what is 
expected for a KA liquid at T = 0.42 as shown in Fig. |3| 

b. 


IV. DISCUSSION 

Changing the stoichiometry in the KA mixture leads 
to important changes in the relaxation dynamics of the 
mixture in the liquid phase. In particular, upon chang¬ 
ing composition towards the 2:1 mixture, the system be¬ 
comes more fragile. Richert and Angell^ suggested that 
more fragile glassformers should be expected to exhibit 
a greater change in structure upon supercooling, as evi¬ 
denced by the jump in specific heat capacity observed in 
experiments on glassforming liquids upon falling out of 
equilibrium at Tg. While there are exceptions to this 
idea such as kinetically constrained models^, spheres 
in higher dimension^ and certain molecular glassform- 
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FIG. 7. Growth rate of the 2:1 KA crystal as a function of 
the temperature, measured by linear fits to the growth curves 
in Fig. [6] The dashed line is a linear fit to these points. The 
crossing point between the linear fit and the zero-line delivers 
an estimate for the melting temperature Tmeit ~ 0.447(2). 


ersSs— , many materials do follow this trend. In addition, 
recent work suggests a link between the size of correlated 
regions and fragility^, supporting the idea that glass- 
formers in which the structure is more developed (involv¬ 
ing more cooperation in the relaxation) are more fragile. 
In previous work on the Kob Andersen mixture and other 
model systems, the degree of structural change is corre¬ 
lated with fragilit y^^i^^i^^ which is corroborated by recent 
work investigating the extent to which the structure pre¬ 
dicts the dynamic a^^d° . 

With the above discussion in mind, our results which 
show the population of LFS decreasing upon changing the 
composition to the more fragile 2:1 mixture for a given 
supercooling [Fig. [21[b) inset], seem odd at first sight. 
We cannot rule out the possibility of the formation of 
some other LFS, but we note the bicapped square an¬ 
tiprism would reasonably satisfy the criteria for a locally 
favoured structure given that the AI 2 CU crystal is stable. 
Furthermore, the transition between the liquid and the 
crystal appears quite strongly first order or at least there 
must be a strong free energy barrier to nucleation as we 
never observe it. 

Intuitively, one might expect the freezing to be weakly 
first order since the LFS matches the crystal structure. 
An example of this is in the 2d hard disc system which 
forms a hexatic crystal in a weakly first order process 
upon increasing the volume fractio n^^d^ . However, we 
find no evidence of homogenous nucleation in our systems 
and conclude that they must be more strongly first order 
than considerations from d =2 would suggest. 

We believe that the explanation lies in the competition 
between different varieties of LFS, each of which con¬ 
tributes differently to the stability and lifetimes of the 
overall LFS population. As schematically illustrated in 
Fig. [H the three main types of bicapped structures differ 
in the role played by A (large) and B (small) particles 


along the principal rotational axis of symmetry of the 
bicapped antriprism. The identity of the capping par¬ 
ticles is crucial, the structures with more of the larger 
A particles have a larger radius of gyration. The radii 
of gyration of AgBa and A 9 B 2 clusters are 4% and 8 % 
smaller than the AiqBi. Therefore, the volume per clus¬ 
ter (and so the volume accessible via rotations around the 
principal axis of symmetry) is also larger for the AiqBi 
and A 9 B 2 clusters than the AgBg cluster. The increase 
in accessible volume and the abundance of A particles 
with respect to B particles in all the considered mix¬ 
tures, make the two energetically more costly structures 
(Fig.[3lD) more favourable from an entropic point of view. 
This behaviour is similar to that observed for a different 
set of competing clusters in the Wahnstrom mixture, an 
alternative model of colloidal suspensions^^. Together, 
the entropic contributions and the ground state energies 
of the clusters, make nucleation hardly accessible on one 
hand and on the other, seeded crystallisation possible. 

We propose, that the AgBa LFS form a subset of LFS 
and that these may play a key role in slowing down the 
dynamics. FigureUKb) inset implies that the 2:1 mixture 
exhibits a stronger dynamic response to the change in 
structure and this composition has a higher population 
of AgBg LFS. We therefore speculate that to understand 
the role of LFS, one should also consider not only its 
geometry but also its composition. If we assume that 
the AgBg is the key LFS, then the original hypothesis is 
satisfied in the sense that the “domains” of crystallisable 
LFS are larger in the 2:1 mixture than in the 4:1. Despite 
the fact that nucleation is not accessible in our compu¬ 
tational time-scales, we expect that were the 2:1 mixture 
to be more deeply supercooled, then the original hypoth¬ 
esis might be satisfied, that at very deep supercooling, 
large domains of AgBg LFS would form. Whether such 
a deeply supercooled state point is sufficiently stable to 
crystallisation to observe these domains of LFS is open 
to question but in any case lies outwit our computational 
resources. 

Any state with sizeable LFS domains would be highly 
supersaturated. Given our observation of crystal growth 
at moderate supercooling, it thus seems likely that were 
any large domains to form the system would then crys¬ 
tallise. We therefore conclude that the relief of frustra¬ 
tion in the KA model by changing its composition may 
be affected by the need to supercool very deeply to gen¬ 
erate large domains of “crystalline” LFS so the system 
would transition to a crystal. In this respect, it may 
be that the KA model satisfies the concept of geometric 
frustration via the change of composition rather than by 
curving space as is usually invoked. 

Frustration has been successfully tuned in some d=2 
systems by curving space^Zr— and other methods^iiS. In 
d =2 the transition to crystallisation is much more weakly 
first order than seems to be the case here, so the LFS pop¬ 
ulation in liquid can be much higher. Conversely, in d=3, 
control of frustration with an LFS compatible with the 
crystal, has been suggested in hard spheres and related 





7 


systems^^— , but the prevalence of fivefold symmetric 
structures may complicate the situatio n^^d^ . 

Finally we note that in certain lattice models more 
frustration against certain types of crystallisation can 
lead to a more complex crystal structure and more weakly 
first order freezing transition. An intriguing avenue for 
future research is how such observations translate in 
atomistic models such as that considered herein. 


V. CONCLUSION 

In summary we have investigated the nature of geomet¬ 
ric frustration in the Kob-Andersen model. Inspired by 
the observation that the 2:1 stoichiometry is compatible 
with the AICI 2 crystal whose unit cell is consistent with 
locally favoured structure of the KA model (the bicapped 
square antiprism), we have explored the relief of geomet¬ 
ric frustration via changing composition rather than by 
curving space. At first sight our results are surprising, 
rather than a strong increase in LFS population upon 
moving to the 2:1 composition, we find a small decrease. 
This is made more surprising by the observation that 
the 2:1 mixture is considerably more fragile than the 4:1 
mixture, which is sometimes correlated with a relation¬ 
ship between structure and dynamics in the supercooled 
liquid. 

We rationalise our observations by considering that the 
LFS themselves can play a role in the behaviour of the 
system. In particular, we have shown that for a given 
LFS population, those LFS whose composition is com¬ 
patible with that of the AI 2 CU crystal, i.e. AgBa are cor¬ 
related with slower dynamics. The population of these 
crystal-compatible LFS indeed increases markedly when 
the composition is changed to the 2:1 mixture, but re¬ 
mains a fraction of the total LFS population. Thus we 
do not find a large increase in LFS domain size upon 
changing the composition. Given that the LFS popu¬ 
lation overall increases upon supercooling, we speculate 
that a more deeply supercooled KA mixture than our 
simulations permit, may yet behave as expected. 
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